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TRANSVERSE  SHEAR  STRESSES  IN  CIRCULARLY 


CURVED  BARS 

1  2 

By  Chlrasak  Thasanatorn  and  Walter  D.  Pilkey  ,  M.  ASCE 
Introduction 

The  exact  stress  determination  of  a  circularly  curved  bar  subjected 
to  pure  bending  and  shear  forces  in  the  plane  of  curvature  requires  the 
solution  of  a  boundary-value  problem  of  three-dimensional  elasticity.  This 
presents  formidable  complications  because  of  the  difficulty  of  satisfying 
the  boundary  conditions.  Few  elasticity  problems  for  curved  bars  have 
been  solved.  Analytical  solutions  of  a  curved  bar  with  circular  cross 
sections  were  developed  by  Golovin,  as  presented  in  Timoshenko  (2).  Also 
given  in  Reference  2  are  plane  elasticity  solutions  for  curved  bars  of 
small  width.  A  so-called  strength-of-material  theory  of  curved  bars  was 
developed  by  Winkler  using  assumptions  analogous  to  the  Bemoulli-Euler 
hypothesis  of  straight  beam  theory.  In  this  case,  the  equilibrium 
conditions  are  satisfied  in  an  average  sense. 

This  paper  presents  a  method  for  the  stress  analysis  of  a  circularly 
curved  bar  subjected  to  a  transverse  end  force  in  the  plane  of  the  bar  curva¬ 
ture.  This  formulation  follows  the  semi-inverse  method  of  Saint  Venant, 
using  a  stress  function.  The  assumption  that  stresses  on  any  cross 
section  depend  upon  forces  at  that  particular  cross  section  permits  the 
three-dimensional  elasticity  problem  to  be  reduced  to  determining  the 
stresses  on  the  cross  section  of  the  curved  bar.  The  stress  solutions  obtained 
from  this  structural  model  may  be  extended  to  more  general  cases  of 
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loading  by  appealing  to  Saint  Venant's  principle,  provided  that  the  cross 
section  of  interest  is  sufficiently  far  from  any  points  of  rapid 
variation  in  the  internal  forces.  Finally, a  finite  element  method,  which 
is  suitable  for  a  curved  bar  with  complicated  cross-sectional  geometry,  is 
used  to  determine  the  stresses  on  the  cross  section. 

Torsional  stresses  for  circularly  curved  bars  are  treated  in 
Reference  1. 
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ELASTICITY  FORMULATION 

Consider  a  planar  curved  bar  of  initial  curvature  R,  with  one  end 
fixed  and  the  other  end  free, as  shown  in  Fig.  1.  The  bar  is  subjected  to  a 
transverse  radial  force  passing  through  the  shear  center  at  the  free  end. 
Assume:  (1)  the  material  is  linearly  elastic,  continuous,  homogeneous,  and 

isotropic;  (2)  the  curved  bar  is  planar,  lying  in.  the  x-z  plane  with 
the  longitudinal  axis  coinciding  with  the  geometric  centroid  of  the  cress 
section;  (3)  the  curved  bar  can  have  any  arbitrary  cross-sectional  configur¬ 
ation  but  must  be  prismatic;  (4)  the  body  forces  are  small  compared  to  the 
magnitudes  of  the  stresses  and  can  be  neglected;  (5)  the  stresses  at  any 
cross  section  are  assumed  to  depend  upon  only  the  internal  forces  at  that 
cross  section;  and  finally  (6)  the  normal  stress  distribution  is  assumed 
to  follow  a  hyperbolic  law  in  the  radial  direction.  The  curved  bar  is 
modeled  with  fixed-free  ends  to  simplify  the  mathematical  formulation.  For 
this  case  the  variations  of  internal  forces  along  the  curved  bar  axis  are 
either  constant  or  are  a  function  of  sine  or  cosine  functions. 

From  the  overall  equilibrium  conditions  of  the  bar  of  Fig.  1,  the 
internal  forces  at  any  cross  section  along  the  bar  axis  are 

P  *  W  sin  0;  *  -TIE.  s±n9;  v  *  W  cos9 

y  (i) 

H  »  M  .  C:  V  ■  0 
x  z  y 

The  twisting  moment  1!^  vanishes  due  to  the  assumption  that  the  force  W 
passes  through  the  shear  center  of  the  cross  section. 

From  assumption  5  (above)  and  Eq.  (1),  the  three-dimensional  elasticity 
problem  may  be  reduced  to  a  two-dimensional  one  by  expressing  the  stresses 
in  the  form 

a  »  a  (y,z)  sin0;  a  -a  (y,z)  sin0;  a  -  a  (v,z)  sin9 
x  x  ^  y  y  y  ’  z  z  '  ’ 

t  *  x  (y, z)  sin6;  t  *  t  (y ,z)  cos9;  t  *  t  (y,z)  coso 
yz  yz  ’  xy  xy  '  xz  xz  J 

where  c.  c  ,  a  ,  t  ,  ~  and  t  are  functions  of  cross-sectional 
x’  y*  z’  yz’  xy  xz 


coordinates  (y,2)  only.  Assumptions  5  and  6  suggest  trying  to  satisfy 
the  stress  equilibrium  equations  of  elasticity  theory  by  assuming  a  normal 


stress  distribution  in  *■>«  form 
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where  C  ,  C  ,  and  C  are  constants  to  be  determined.  The  stress  expressions 
o  y  z 

_hat  automatically  satisfy  the  three  stress  equilibrium  equations  can  be 
expressed  in  terms  of  the  two  stress  functions  v(y,z)  and  '? (y,z) 
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where  F(z)  and  G(y)  are  defined  along  the  cross-sectional  boundary. 


The  stress  expressions  in  Eqs.  (4)  -  (9)  are  derived  independent  of 
the  displacements  u,  v,  and  w.  The  stress  components,  or  more  specifically 
the  stress  functions  and  V ,  can  not  he  taken  as  arbitrary  functions  of 
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(8,  y,  z)  but  are  subject  to  the  two  stress  compatibility  conditions  to 
ensure  the  existence  of  displacements  u,  v,  and  w  corresponding  to  these 
stresses.  A  set  of  consistent  stress-displacement  relations  can  be 
derived  by  using  the  Hellinger-Reissner  variational  principle,  which  is  a  mixed 
principle.  Then,  two  independent  stress  compatibility  conditions  can  be 
derived  from  the  five  consistant  stress-displacement  equations.  They  are 
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where  r  *  (R  +  z).  Cartesian  and  cylindrical  coordinates  are  used  inter- 

changably  as  convenience  dictates.  Tne  constants  C  ,  C  ,  and  C  in 

o  y  z 

the  assumed  normal  stress  expression  of  Eq.  (2)  are  determined  from  the 

following  equilibrium  conditions  for  any  section  of  the  curved  bar. 
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The  constants  C  ,  C  ,  and  C  are  found  to  be 
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where 
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By  substituting  Eqs.  (4)  -  (9)  into  (10)  and  (11),  the  two  governing  partial 
differential  equations  for  the  stress  functions  $  and  V  can  be  obtained. 

R2  .  4R2  6R2  R 

(R  +  z)2  ^’zzzz  (R  +  z)3  ’zzz  (R  +  z)4  ^’zz  (F.  +  z)  ’zzzz 


dydz 


2R 


2R 


+  (R  +  z)2  ‘ ’ zzz  (R  +  z)3  ^’zz  +  (R  +  z) 2  d ’ yyyy  (R  +  z)  "’yyyy 


2R 


2R 


l  Rz 


4R2 


(R  4-  z)  ‘ ’ yyzz  +  (R  4-  z) 2  ‘ ’yyz  +  (R  +  z) 2  ^ ’ yyzz  (R  +  z) 3  ^ ’ yy 


2R2 

V  v  . - \~2  4>, 


2R 


4R2 


|  (R  4-  z)z  v,yy z?  (R  4-  z)  ‘’yyzz  (R  +  z)3  ^ ’ yyz  (R  +  z)1*  ^’yy 


6R2 


yvz 


2R" 


2R 


2R 


2R 


(R  4-  z)2  ’yyz  (R.  4-  z)3  ’’yy  (R  +  z)  “’yyzz  (R  4-  z)z  *’yyz 

4R2 


2R2  _ 

(R  +  z)2  ’’’yyzz  (R 


>r2  a  I  ,  fi  er; 

4-  z)3  ?’yyz  j  L2  (R  + 


2  2  (F(z)  -  y) 


z)  ~  dz 


(Co  *  2f7m»],2I  +  I  [f^y  +  Cz  G(y) 

(c0  +  2V, 


d-G(y) 
x  — — + 
dyz 


r 1  er2  fc 

/  g(yj  i 

\  .  C  /  G2(y) 

2P-C(y) 

|_2  (R  +  z) 

a  \  (R  4-  z) 

)  Cz  \(R  +  z) 

(R  +  t)  R  j 

’  zz 

+  ff?iT3  (Co  +  2Cyy)  +  TF+"z)3  (Co  +  V  •  c*a!  ■  0 


(16) 


7 


3R2 


(R  +  z)2  ?’zzz  (R  +  z)3  ^’zz  (R  +  z)  ‘’zzz  (R  +  z)2  ‘  ’zz 


V 


+  7U2R  T? 


+  v ! 


1  .  i  ,  3R2  R  „  2R  | 

j  (R  +  z) 2  ’vyz  (R+z)3  ’yy  (R+z)  ‘ ’yyz  (R+z) 2  *yy  j 


+  (1  +  v)  !  - 


V 


,  R2  , 

+  T2  <P  » , 


|  (R  +  z)  ''yyz  (R  +  z)2  ^’yyz  (R  +  z)  ^’yyyy  (R  +  z)  ^ ’ yyzz 

Sy*  |+  <R  +  z)[t  tfvzT3  tt1  (f(2)  -  y)<co +  v(*>)]-. 


R2 


(R 


_  (R  +  _jl  EvRj  ,  fc  ,  G(y)  _  n  ,  r  ,  G2 (v) 
(R  Z)[2  (R  +  z)2  LCo  (R  +  z)  1  Cz((R  +  z) 


2RG(v) 
(R  +  z) 


+  R  -  z) 


J]- 


(17) 


2EvR2  / 

(R+z)2  \ 

'T  aTTIT  Khco  +  c2(«y>  -  2R>]-5 


[°o +  y +  h<z  - R)) 


+  a  +  ">  (co  +  2V> 


The  traction-free  boundary  conditions  on  the  lateral  surface  of  a  curved 


bar  yield  the  values 
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for  the  stress  functions  along  the  cross-sectional  boundaries. 
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FINITE  ELEMENT  ANALYSIS 


The  governing  equations  for  the  stress  functions  $  and  V  as  given  in 
Equations  (16)  and  (17)  are  too  complex  to  be  suitable  for  closed  form 
solutions.  A  variatiui^lly  based  finite  element  procedure  will  be  applied 
to  obtain  the  solutions.  The  problem  of  integrating  the  partial  differen¬ 
tial  equations  (16)  and  (17)  subject  to  boundary  conditions  (18)  may  be 
transformed  into  the  problem  of  finding  the  unknown  stress  functions  <$> 
and  V  which  make  stationary  a  functional  subject  to  the  same  boundary 
conditions.  The  functional  may  be  expressed  in  matrix  form  as 
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The  procedure  will  be  to  idealize  the  cross-sectional  domain  by 
an  assemblage  of  triangular  elements.  The  functional  H  may  be  approximated 
as  the  sum  of  functionals  evaluated  for  each  element.  Mathematically, 
the  functional  may  be  written  as 
M  (e) 

II  -  E  -  H  (24) 

i-1 


where  i  indicates  the  element  number  and  M  is  the  total  number  of  elements. 
At  node  j  the  stress  ^unctions  and  their  first  derivatives  (slopes)  with 
respect  to  y  and  z  are  chosen  as  the  unknown  nodal  parameters  and  are  ex¬ 
pressed  in  matrix  form  as 


The  stress  function  behavior  within  each  element  may  be  written  in 
terms  of  their  nodal  parameters 
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where  [N]  is  a  shape  function  matrix  and  is  a  vector  of  the  unknown 

nodal  parameters.  The  vector  of  derivatives  in  Eq.  (20)  may  be  obtained 

by  differentiating  Eq.  (26)  with  respect  to  y  and  z.  Thus 
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where  elements  of  matrix  [B]  are  the  second  derivatives  of  the  shape  func¬ 
tions  with  respect  to  y  and  z  as  defined  in  Eq.  (20).  The  discretized 
functional  for  element  i  is  obtained  in  terms  of  the  unknown  nodal  para¬ 
meters  by  substituting  Eq.  (27)  into  (19) 
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where  integration  is  performed  over  the  element  domain  A^.  Finally  the 
total  functional  is  obtained  by  collecting  the  contributions  from  all 
discretized  functionals  within  the  cross-sectional  domain. 
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The  stationary  solution  of  the  functional  I!  can  be  derived  from 
the  condition  that  the  variation  of  the  functional  J!<y,z,i»^)  with  respect 
to  the  unknown  nodal  parameters  vanishes .  That  is 
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where  N  is  the  number  of  nodes  on  the  cross  section. 
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NUMERICAL  RESULTS 

The  formulation  will  be  applied  to  two  curved  bar  cross  sections. 
First,  a  homogeneous  thin  curved  bar  of  rectangular  cross  section  will 
be  used  so  that  the  calculated  stresses  can  be  compared  with  an  analytical 
solution  available  from  plane  elasticity  theory.  The  second  example,  is 
chosen  to  be  a  curved  box  section. 

Homogeneous  Thin  Curved  Bar  of  Rectangular  Cross  Section.  The 
thin  rectangular  curved  bar  of  Fig.  (2a)  has  a  1  in.  (0.0254  m)  width,  10 
inch  (0.254  n)  depth,  and  a  radius  of  curvature  of  20  inches  (0.508  m) .  The 
cross  section  is  subjected  to  an  internal  shear  force  in  the  radial  direction 
of  10,000  lb  (44.48  kN)  at  its  centroid.  The  discretized  model  of  the  cross 
section  is  shown  in  Fig.  (2b)  and  partial  computer  results  are  presented  in 
Fig.  (3).  The  stress  results  obtained  by  the  present  formulation  are 
compared  in  Fig.  (4)  with  those  obtained  from  the  plane  elasticity  theory. 

The  direct  shear  stresses  based  on  the  present  formulation  are  the  average 
values  across  the  bar  width. 

Homogeneous  Curved  Box  Girder  Subjected  to  Direct  Shear  Force. 

The  second  problem  demonstrates  the  applicability  of  the  method  of  analysis 

to  a  more  complicated  cross-sectional  configuration.  Choose  the  curved 

bar  section  of  Fig.  (5).  'it  is  made  of  steel  with  a  modulus  of  elasticity 

6  2 

and  Poisson's  ratio  of  30  x  10  psi  (206.85  GN/m4 )  and  0.3  respectively. 

The  section  is  subjected  to  an  internal  shear  force  of  100,000  lb  (444.8  kN) 
in  the  z  direction.  The  finite  element  idealization  of  the  cross  section 
is  shown  in  Fig.  (5).  The  shear  stresses  throughout  the  walls  were  computed. 
The  average  shear  stresses  across  the  wall  thickness  are  illustrated  in 
Fig.  (6). 


CONCLUSION 


The  paper  presents  a  formulation  for  the  direct  (transverse)  shear 
stress  analysis  of  a  circularly  curved  structural  member.  The  formulation 
is  based  on  the  three-dimensional  theory  of  elasticity.  The  assumption 
"stresses  at  any  cross  section  depend  only  on  the  internal  forces  at  that 
particular  cross  section"  enables  the  three  dimensional  elasticity  theory 
to  be  reduced  to  the  determination  of  the  stresses  at  a  particular  (two- 
dimensional)  cross  section.  Although  the  solution  is  based  on  a  cantilevered 
bar  with  a  particular  applied  loading,  this  formulation  can  be  used  to  compute 
the  transverse  shear  stresses  for  a  circularly  curved  bar  with  any  boundary 
conditions  and  arbitrary  loading.  Thus,  the  formulation  requires  only  that 
the  internal  forces  be  known  on  the  cross  section  where  the  shear  stresses 
are  sought. 
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Fig.  3  Direct  Shear  Stress  Distribution  (in  psi) 
on  the  Thin  Curved  Bar  Section  of  Fig.  2 

Note:  1  in.  *  0.0254  m:  1  psi  ■  6.895  kN/cT 
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Fig.  4a  Comparison  of  the  Shear  Stresses  (in  psi)  along 
the  Center  Line  of  a  Thin  Curved  Bar  of  Fig.  2 
Obtained  by  the  Present  Formulation  with  those 
of  Two-Dimensional  Plane  Stress  Theory 


Fig.  4b  Direct  Shear  Stress  Distributions  Along  the  Center 
Line  of  the  Bar  of  Fin.  2 
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